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ABSTRACT 

Star clusters that form in nuclear rings appear to be at slightly larger radii than the gas. We argue 
that the star clusters move out from the gas in which they are formed because of satellite-disk tidal 
interactions. In calculating the dynamics of this star cluster and gas ring system, we include the effects 
of dynamical friction of the background stars in the host galaxy on the star cluster, and inflowing gas 
along the bar onto the nuclear ring at the two contact points. We show that the final separation is 
of the order of the Hill radius of the nuclear ring, which is typically 20 — 30% of its radius. Massive 
star clusters can reach half of this separation very quickly and produce a factor of a few enhancement 
in the gas surface density. If this leads to star formation in addition to the (ongoing) formation of 
star clusters near the contact points, a possible (initial) azimuthal age gradient may become diluted 
or even disappear. Finally, if the star cluster are massive and/or numerous enough, we expect the 
nuclear ring to migrate inward, away from the (possibly) associated (inner) Lindblad resonance. We 
discuss how these predictions may be tested observationally. 

Subject headings: accretion, accretion disks — stellar dynamics — stars: formation — galaxies: bulges 
— galaxies: star clusters — galaxies: nuclei 



1. INTRODUCTION 

Nuclear rings in barred spiral galaxies are regions of 
large gas surface densities and high star formation. Their 
existence is intimately related to the stellar bars with 
which they are associated. The stellar bar drives the 
interstellar medium (ISM) into the nuclear ring. How 
the radius of the ring is determined is unclear. The 
radius may be set by the location of the inner Lind- 
blad r esonance (ILR) of the stellar bar (see review by 
iButa fc Combes 1996, and references therein). The rea- 
soning behind this idea is that the sign of the torque 
changes on either side of the Lindblad resonance. Alter- 
natively, the radius may be at the location where the non- 
axisymmetric bar weak ens, which is typic ally at ~ 10% of 
the size of the bar (Shlo sman et al.l[l990( ). Finally, it may 
also be determined by where the transition from phase 
space domination of Xi orbits (parallel to the major axis 
of the large-scale bar) to the phase space d omination of 
X2 or bits (perpendicular to the bar) occurs (| Regan et al.l 
[1991 . 

Regardless of how they are formed they are interesting 
as a class because they are perhaps the largest popu- 
lation of nearby star bursting regions. They form star 
clusters in prodigious amounts and are the only environ- 
ment where super star clusters (SSCs) might be found in 
abundance in normal galaxies (jMaoz et al.ll2001l ). These 
SSCs can have masses in excess of 10^ Mq and may be 
the progenitors of globular clusters if they can survive 
for 10^° years. 

In some of these systems the location of these star clus- 
ters is curious. They appear to be at larger radii than 
the gas ring from which they are presumably formed. 
For instance, the morphology of NGC 1512 clearly show 
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the star cluster(s) exterior (radially) to the gas as op- 
posed to being embedded in the gas (jMaoz et al.1 12001). 
Taking another example, in NGC 4314 the star clus- 
ters are at larger rad ii than t he gas in the nuclear ring 
(jBenedict et alj|2002t) . Finally. [Martini et all (|2003[) find 
that in their sample of 123 galaxies, eight have strong nu- 
clear rings and in each of the eight, star formation occurs 
outside of the dust ring. Two of their other galaxies have 
dust lanes immediately interior to what appears to be an 
older stellar ring. 

There are two possible explanations for why these star 
clusters are external to the gas. First, while the gas can 
lose angular momentum due to dissipation, the stars and 
star cli isters are coUisionless a nd hence remain on larger 
orbits ([Regan fc Teubenll2003f) . Second, the star clusters 
migrated outward after their formation. We explore the 
latter possibility in this paper. 

We argue that tidal interactions between the star clus- 
ter and the nuclear ring from which it was formed leads 
to their separation. We apply the physics of "gas shep- 
herding", which was studied by one of us ([Chang|[2008L 
hereafter COS), to this cluster- ring system. We derive a 
natural length scale for this separation, which we call the 
"Hill radius" of the ring. At the same time, we show that 
the outer edge of the ring can experience a surface den- 
sity enhancement of a few, possibly leading to additional 
formation of stars and star clusters along the edge of the 
ring. Finally, we find that if the star clusters are mas- 
sive and/or numerous enough the cluster-ring systems as 
a whole is expected to migrate inward. We show that 
these predictions can be tested observationally. 

We present our model in § [5| After summarizing the 
observed properties of nuclear rings in § 12.11 we give the 
basic picture and introduce the Hill radius of the ring 
followed by the basic equations of our model in 



2.31 The solutions for various scenarios are outlined in 
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§|3l A few observationally testable predictions that arises 
from our model are given and compared to observations 
in § m We discuss some aspects of our model in § [5] and 
summarize our conclusions in § [6l 

2. TIDAL INTERACTIONS BETWEEN STAR CLUSTERS 
AND NUCLEAR RINGS 

We first summarize the basic properties of observed 
nuclear rings, which are used to place the model predic- 
tions in context. We then present the basic picture of our 
model, followed by the basic set of equations that govern 
our model. 

2.1. Observed properties 

For a sample of 22 nuclear rings. Mazzu ca et al.l ()2008f ) 
find that the rings are in the same plane as the disk of 
their host galaxy and that they are intrinsically circular. 
Typically, the radius of the ring is around ~ 0.5 kpc, but 
it can be a factor three smaller or larger. With a circular 
velocity at the ring radius of around ~ 150kms~^, the 
(spherically) enclosed mass is about Monc ~ 2.5x 10® M©. 

The gaseous mass in the ring may be estimated from 
high spatial resolution CO measurements, using the 
'standard' conversion to H2 mass and adding 36% for 
He. In this way, we obtain fo r NGC5457 NGC 3351, 
NGC 6951. NGC 3504 (KennevlTal] [1991 [l993l) and 
NGC4314 (jBenedict et al.lll996[) on average a mass in- 
side the ring of M^ ~ 5x 10^ Mq. Although the ring mass 
may be as much as a factor ten smaller or two larger, 
it correlates with the ring radius, so that the ring-to- 
enclosed mass ratio g^ = M^/Menc stays within a narrow 
range of qd = 0.01 — 0.03. While this is the molecular 
gas inside the ring, part of the gas is being ionized, due 
to photons from newly formed massive stars or due to 
shocks from stellar winds or supernovae explosions, and 
at the inner edge of the ring possibly even due to an 
active galactic nucleus (AGN). However, the ionized gas 
mass estimated from the Ha (plus [N II]) emission is gen- 
erally two ord ers of magnitude low er than the molecular 
gas mass (e.g. IPla nesas et a"Llll997[ ). 

The nuclear rings are easily observable due to the emis- 
sion from the massive stars inside the young star clusters 
in the ring (which in the ultraviolet even dominates the 
emission from a possible AGN) . The combination of high- 
spatial resolution (HST) images and multi-wavelength 
coverage, makes it possible to isolate the individual clus- 
ters and after correction for (dust) extinction to estimate 
th eir mass and (re l ative) age. Below we use such a study 
bv lBenedict et al.l (|2002D to investigate in more detail the 
nuclear ring in NGC 4314. Overall, the cluster masses 
follow a power-law distribution with index —2, similar 
to young s tar clusters in mer ging galaxies such as the 
Antennae (jZhang fc Falll [l999l ) . but typically extending 
to a smaller upper limit in mass (few times 10^ Mq), 
i.e., fewer of the so-called "super star clusters" (but see 
iMaoz et al]|200lD . 

On the other hand, these young star clusters are most 
likely just the latest s tar formation in a n ongoing suc- 
cession of bursts (e.g. lAllard et al.l 1200(f ). Most of the 
lower-mass clusters dissolve, either early-on (< 50Myr) 
due to violent relaxation after expelling the residual gas 
('infant mort ality') depending on th e star formation ef- 
iiciency (e.g. iParmentier et al]|2008[ ). or later-on due to 
internal relaxation ('evaporation') and external tidal ef- 



fects fe.g. [McLaughlin fc Falill2007[ ). The resulting evo- 
lution towards a bell-shaped distribution implies a rel- 
ative increase in more massive star clusters over time. 
Unfortunately, once the most massive stars have died, 
the luminosity of star clusters very rapidly fades, so that 
it becomes more difficult to distinguish them from the 
old stellar population of their host galaxy. Nevertheless, 
based on velocit y dispersion measu rements from stellar 
absorption lines. iHagele et all (|2007f ) estimate dynamical 
masses up to almost ^ 10^ Mq for individual star clus- 
ters in the nuclear ring of NGC 3351. Henceforth, the 
typical ratio of cluster-to-ring mass ratios q' = M^/Md 
depends on the (average) age of the star clusters. For 
the young, luminous star clusters we adopt the range 
q' = 0.001 — 0.01, whereas for the old(er), faint(er) star 
cluster q' = 0.01 — 0.1 is probably more appropriate. Be- 
low we investigate the two cases q' = 0.1 and q' ~ 0.01, 
corresponding to M^ ~ 5 x 10*^ M© and M^ ~ 5 x 10^ Mq 
in case of the above average gas mass in nuclear rings of 
Md ~ 5x10^ Mq . The effect of a less massive star cluster 
with q' < 0.01 is similar to that of q' = 0.01. 

The time for a star cluster to cross the nuclear ring is 
given by icross = 2n/D,, where 51 — ^GMcnc/f^ is the 
orbital frequency and G is Newton's constant of gravity. 
Given a typical circular velocity of ~ 150 km s^^ at the 
observed mean nuclear ring radius of 0.5 kpc, the cross- 
ing time is around icross ~ 20Myr. A star cluster is 
subject to dynamical friction with the bulge stars on a 
(minimum) timescale^ idf = -/Wenc/(^-Ws In A), with In A 
the Coulomb logarithm. Using tcross = ^df 'z'o'd 27rln A 
and substituting qd — 0.02 from above, we find that the 
dynamical friction time in units of the crossing time is 
approximately given by tdf/^cross — l/*?'- This means 
hundreds to tens of orbits for star clusters with masses 
that are 1% {q' = 0.01) to 10% (g' = 0.1) of the nuclear 
ring mass, or, with icross = 20Myr, a dynamical friction 
time tdf from about 2 Gyr to 200 Myr for a non-rotating 
bulge, which we will assume for the remainder of this 
work. 

2.2. Basic Picture 

We illustrate our basic model in a cartoon in Figure 
[TJ Gas fiowing in a bar potential along the Xi orbits 
transitions to X2 orbits making up the nuclear ring. The 
X2 orbits are approximately circular and so we will ap- 
proximate them as circles. From the gas in the nuclear 
ring stars are formed in star clusters. We assume that 
the star clusters are formed preferentially near the outer 
edge. These newly formed and forming star clusters are 
subject to dynamical friction with the background bulge 
stars and tidal interactions with the ring. However, the 
initial tidal interactions are strong compared to dynami- 
cal friction and thus these star clusters will initially move 
away from the ring. 

Indeed, in observations of nuclear rings young star 
clusters often appear to be at larger radii than the 
gas ring from which they are presumably formed. The 

^ The dynamical friction time is likely to be longer for a partially 
rotationally supported bulge as opposed to a non-rotating bulge. 

■* The background bulge stars may also be halo stars, and the 
bulge/halo is assumed not to be supported by rotation. For the 
purposes of this work, we will idealize the bulge/halo as a n isother- 
mal background distribution, which we argue in § 15.11 is a valid 
approximation. 
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color comp osite images of NGC 1 512 (|Maoz et al.ll2001[ ). 
NGC4314 (jBenedict et all 12001) . and NGC 7742 (Hub- 
ble Heritage^), show young star cluster as bright spots on 
the outside of gas traced by the dust. As pointed out by 
the referee, from these images it is also clear that, even 
though the nuclear rings are as we assume approximately 
circular with star clusters at the outer edge, the observa- 
tions can be quite difficult and complex. First, when re- 
solved in the optical the rings can look like tightly wound 
spirals. This might be just appearance because the bar 
dust lanes connecting with the ring are also partly ob- 
scuring the ring in optical broad-band images, whereas in 
e.g. narrow-band imaging around Ha the ri ng of young 
star clusters seem to be complete (see also iMaoz etaLl 
l2001f) . Nevertheless, our basic picture described above 
for a circular ring should still provide a fair description 
in case of a tightly wound spiral. Second, very obscured 
star clusters could remain hidden below the gas ring even 
in the near infrared images taken, whereas even longer 
wavelength observations (yet) lack the high spatial reso- 
lution required to resolve the star clusters. Even so, the 
clearing of gas and dust from the newly formed star clus- 
ter is very fast and efficient, as even the you ngest visible 
star clusters are only mildly reddened (e.g. iMaoz et al.l 
|2001|) . Finally, some (sectors of) star clusters seem to 
appear more on the inside of the nuclear ring, possibly 
along an inward extension of the gas inspiraling onto the 
ring along the large-scale bar. As we discuss briefly in 
§ 15.31 star clusters may form within the ring, in par- 
ticular at the contact points where the gas flows into 
the ring. Moreover, the interaction of non-axisymmetric 
structures inside the ring (such as nuclear bars and spi- 
rals) with the interstellar medium may also induce the 
formation of star clusters. However, even though such 
additional mechanism might be operating and the ob- 
servations should be interpret with care, they in general 
support our basic picture of a circular nuclear ring with 
young stars migrating outward. 

Since the basic physics mirror the "gas shepherding" 
scenario studied by COS, we adopt the same notation. 
In this case the satellite is a star cluster and the disk is 
a ring, but since only the (local) geometry at the edge 
matters the physics remains the same. There are however 
two important differences with respect to COS. First, 
the perturbing star cluster's initial radial position begins 
very close to the ring as opposed to a satellite reaching 
the disk from very far away in COS. Due to ring tidal 
interactions, the star cluster and ring will separate to 
establish an initial separation, Ar, which we show below 
is of order the "Hill radius" of the ring. We aim to study 
the manner by which this initial separation is achieved. 
Second, gas continues to flow in from the Xi orbits so 
that at the transition between the Xi to X2 orbits, gas 
is continuously being added to the system. Therefore, a 
source of mass exists in this system, which changes the 
dynamics. 

For now we ignore the source of mass to concentrate 
just on the tidal interaction between the ring and star 
cluster and the dynamical friction experienced by the 
star cluster. We begin by showing that this combined 
action of dynamical friction and ring tides results in a 
length scale in the problem, which we call the Hill radius 

5 http://heritage.stsci.edu/1998/28/ 



of the ring. T he torque due to dynamical friction on the 
star c luster is (|Chandrasekhail[l943l : iBinnev fc Tremaini 
flQSl 



Tdf 



hr^s 



(1) 



where Mg is the mass of the star cluster, rg is its dis- 
tance from the center of the galaxy, Og is its orbital 
frequency, and, as before, t^i is the timescale for dy- 
namical friction. The torque due to the excitation of 
spira l den sity waves in a disk i s (Goldreich fc Trcmai^ 
19801: [Lm fc' Papaloizou 19S6t IWard fc HouriganI [19891 : 



Artvmowicz..l99a: Ward.l997f) 
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where Md is the mass of the ring with radius r^. The 
radial separation between the star cluster and the ring 
Ar = Ts — rd is assumed to be small compared to the 
distance to the center, i.e., Ar <C r^ < r^. While dy- 
namical friction torques down the star cluster, the ring 
torques up the star cluster. Balancing these two torques, 
i.e., Tdf ^ Tdisk, gives 



Ar 



/ Md 



\M, 



1/3 



rH,d, 



(3) 



where the latter defines the Hill radius of the ring. 

Given the observed range in ring-to-enclosed mass of 
qd = 0.01-0.03 (§[211]), the HiU radius is around 20-30% 
of the radius r^ of the star cluster. Since r^ ^ r^ 3> 
rn.d, the observe d range in nuclear r ing radii of r^ = 
0.2- 1.7kpc fe.g. lMazzuca et al.l[2008[ ) results in tha = 
50 — 500 pc. We expect the radial separation between the 
star cluster and nuclear ring to be of order of this Hill 
radius. However, the exact amount and nature by which 
this separation will be achieved depends strongly on the 
strength of the tidal torque between the cluster and the 
ring, which we intend to quantify in this paper. 

The second process which we wish to study is the effect 
of an additional source of mass on the dynamics of the 
cluster-ring system. The gas merges onto the nuclear ring 
at the two contact points where the Xi orbits transition 
to X2 orbits. We assume that the region where the gas is 
added to the nuclear ring is radially narrow compared to 
the radius of the nuclear ring. For a mass inflow rate of 
M, the amount of angular momentum added to a system 
per unit time is 

L^MVlr^. (4) 

However, not all of this angular momentum can be im- 
mediately applied to affecting the qualitative dynamics 
of the system. By qualitative, we mean that the mass in- 
flow rate, which acts as a source of angular momentum, 
changes the behavior of the cluster-ring system. For in- 
stance, the source of angular momentum from mass in- 
flow balances the sink of angular momentum from dy- 
namical friction on the star cluster. To estimate what 
mass inflow rate would be needed to effect such a change, 
we suppose this mass is added just inside of the satellite 
radial position, i.e., r ~ 7-5. The mass will then be pushed 
to the edge of the ring, rd^ due to tidal interactions. The 
change in radial position will be of order the Hill radius 
of the disk, tha- Hence, for such a parcel of gas of mass 



van de Ven & Chang 

Gas Flow along Xi orbits 



Forming star clusters 




Nuclear Ring in X2 orbits 



Gas Flow along Xi orbits 



Ar ~ rH,d 



Fig. 1. — Cartoon of our model. Gas flows from larger radii (along a X\ orbit parallel to the major axis of the large-scale bar) onto the 
nuclear ring, which is in an approximately circular X2 orbit. Part of this gas collapses to form star clusters near the outer edge of the ring. 
These star clusters, being subject to dynamical friction with the surrounding bulge stars in the host galaxy, tidally interact with the gas in 
the nuclear ring and migrate outward to a radius Ts- At the same time, the star clusters that are sufficiently massive push the outer edge 
of the nuclear ring to a smaller radius r^. The net separation between the star clusters and nuclear ring, Ar = Ts —t^, is of order the "Hill 
radius" of the nuclear ring, r^f d- 

M, the change in angular momentum moving from r^ to 
Td is 

Ai -- M nr Ar -- M Q^rs rH,d- (5) 



For a constant mass inflow rate M, we expect the torque 
from this inflow of gas, acting on the star cluster, to be 

Tinfi = —T— ^ M n^rs THM- (6) 

Setting this inflow torque equal to the dynamical friction 
torque in eq. ((1]), we find the scale for the mass inflow 
rate to be 



M ' 






Mr: 



Md 



1/3 



(7) 



where we have used the definition of the Hill radius of the 
ring given in eq. dH) . Had we instead (naively) set eq. ([4]) 
to be equal to eq. (P), we would have found M ~ Ms/tdf , 
which can be significantly smaller than the more physical 
scale in eq. ([7]). 

As shown in § 12.11 above, idt — Across/?', so that 
M ^ (Md/tcross) I'^/Qd ■ Substituting the average (gas) 
mass inside the ring of Md ^ 'a x 10^ M©, the crossing 
time of icross ^ 20 Myr, and the observed range in ring- 
to-cnclosed mass of qd = 0.01 — 0.03, the typical scale 
for the mass inflow rate becomes M ~ lOg'^ M© yr~^. 
For star clusters with masses that are 1% (g' = 0.01) 
to 10% {q' = 0.1) of the nuclear ring mass, this im- 
plies a scale M from about 0.001 to O.lMoyr^^ To 
balance the combined dynamical friction from the many 
star clusters that are formed, the total mass inflow rate 
might easily have to be as high as ^ lM0yr~^, which 



is the typical star f ormation rate in nuclear rings (e.g. 
iMazzuca et al.ll2008[ ). For mass inflow rates signiflcantly 
below this scale, we expect the effect of this source of 
(gas) mass to be insignificant. For mass inflow rates at 
or above this scale, the effect of dynamical friction to 
push the cluster-ring system as a whole to smaller radii 
is halted. 



2.3. Basic Equations 

We now present the governing equations for our model. 
We do not present a detailed and complete derivation be- 
cause these equations have been already derived in much 
of the literature (e.g. [L in & Paoaloizou 1979b, a, 198a 
Hourigan & Ward 1984. IWard fc Hour igan 1989; Rafiko3 
2002), with two differences. First, we include the effect of 
dynamical friction of the background bulge stars, which 
is extensively discussed in C08 (to which we refer the in- 
terested reader). Second, we will include a mass source 
term for, which models the inflow of gas mass transition- 
ing from the Xi to X2 orbits. 

The evolution of a viscous gas ring under the influence 
of an external torque is given by equations for continuity 
and angular momentum. The equation of continuity is 



9S 
'dt 



1 d{rVrY.) 
r dr 



= S{r,t), 



(8) 



where r is the radial coordinate, v^ is the radial compo- 
nent of the velocity, E is the surface mass density of the 
gas in the ring, and S{r, t) is a position and time depen- 
dent source term. The angular momentum equation is 
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dr 



dr 
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(9) 
where Tvisc = — 27rr j/E dfl/dr is the viscous torque with 
V the viscosity. Taking into account that Q (through 
A^onc) is not an explicit function of time, we combine 
eq. ([S]) and eq. ([9]) to eliminate Vr, so that we are left 
with 

dt 2'Kr dr dr \ dr dr 

(10) 
The radial motion of the star cluster is governed by the 
competing torques of dynamical friction and ring tides. 
Henceforth, the change of angular momentum of the star 
cluster with time is 



S{r,t). 



KL 



dt 



= T, 



df 



Tc 



disk- 



(11) 



For the (bulge of the) spiral galaxy in which star clus- 
ter and nuclear ring are embedded, we adopt a density 



profile p ex r~ (see § 15.11 for a discussion on the va- 
lidity of such an isothermal profile). The correspond- 
ing mass profile is given by Mcnc = .^cnc.o (''/^o)i with 
-^^cnc,o the mass enclosed within radius tq. Such a 
mass profile has the nice property that the orbital ve- 
locity, Worb = r^r = \/GMcnc,o/ro, is constant, and thus 
9(risrj)/9i — Voihdrs/dt in eq. (fTTjl . With this mass 
profile, the partial derivatives of the viscous and tidal 
torques reduce to (see also COS) 



dr 

9Tdisk 
dr 
viscosity 
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follow COS, 
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(12) 
(13) 



adoptin; 



ig J^ = 

3 (Il986h 



For the 

vo(T,/T,q)^ as suggested by iLin fc Papaloizor 

and vq — MQaQd )^ /'^'Td.oVoih using the standard 

IShakura fc SvunvaevI ([1973)' a-disk with iToomrd (|1964D 

Qo parameter, and, as before, qd = Md/Mcnc,o the ratio 

of the nuclear ring (gas) mass to the enclosed mass. 

Substituting these expressions and rescaling the sur- 
face mass density and time respectively as cr = S/Eq 
and t' = i/idf,o, with tdf,o = Mcnc,o/(^oA^s In A), equa- 
tions pUj) and pT|) become 



da_a' <o52(m3) 



dt' ^ q''^^' 



dr' 



(i'q'qa 



^io d 



so that for 7=1, the inflow torque and dynamical fric- 
tion torque approximately cancel, corresponding to the 
t),€is^ of the mass inflow rate in eq. ([7]) . 

Equations (I14|) and (jlSp are the governing equations 
governing our model. They are the same as respectively 
equations (23) and (24) in COS, except for the extra 
source term at the right-hand-side of eq. (HH). In ad- 
dition, our initial conditions and timescales of interest 
are different from the scenario studied in COS. 

3. RESULTS 

We flrst estimate the rate at which a star cluster moves 
out from the nuclear ring. Next, we numerically solve the 
equations describing our model, without any additional 
source of gas, and in case of a constant mass inflow rate. 

3.1. Estimate of the Rate of Cluster-Ring Separation 

We give an order of magnitude estimate for the rate 
at which a newly formed star cluster will separate itself 
form the nuclear ring. For simplicity, we suppose that 
the gas in the nuclear ring does not react signiflcantly 
to the tidal torque, i.e., S ~ Eg (u ~ 1) and r^ ~ rdfi 
remain approximately constant. In this case, eq. (|14p can 
be neglected, and ignoring the last term in eq. p^ due 
to dynamical friction, it simplifies to 



drs 
dt' 



2(3'qdrd,o 



(rs 



■ dr. 



(17) 



Since rs — r <^ r, the latter integral out to the edge of 
the ring r^^o can be approximated as [rd,o/(^s — i^d.o)]'^/^- 
This leads to an ordinary differential equation in x = 
[rs - rd,o)/rd,Q of the form dx/dt' = (2/3'gd/3) x'^. The 
resulting solution, together with the Hill radius ru.d de- 
fined in eq. ^ and r^ > rdfi ^ 'rH,d, implies that 
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As before, q' = Ms/Md is the ratio of the cluster mass 
to the nuclear ring (gas) mass, a' = a(5o/(21nA) and 
/?' = /3/(27rln A) set the strength of the ring viscous and 
tidal torques relative to the dynamical friction torque, 
and the constant 7' controls the amount of inflow torque 
due to gas coming in along Xi orbits that transitions to 
X2 orbits. The dimensionless function /(r) defines the 
radial extent where this gas is deposited in the nuclear 
ring. We choose the form 



f{r) = 



= / 2S? if ^"i.o - Ar < r < rd,o 



otherwise 



(16) 



Since the left-hand-side as well as 8/3'/3 are of order 
unity, the timescale for a star cluster to separate itself 
from the nuclear ring in which it was formed is of or- 
der tsop ~ {rH,d/rd,o)tdffi- The observed range of qd = 
0.01 — 0.03 implies tscp/^df,o ~ 0-2 — 0.3. Furthermore, 
in units of the crossing time icross — idf,o9'9d27rln A, 
the separation time is about isep/^cross ~ {rH,d/rd,o)/<}' ■ 
This means a few to tens of orbits for cluster-to-ring mass 
ratios from q' = 0.01 to q' = 0.1, or about 500 Myr to 
.ftJ'Myr for a typical crossing time of icross = 20 Myr. 

Since the time t in eq. p8|) goes as a fourth power of 
the cluster-ring separation rg — rd.o, the initial separation 
i^l&^uch faster than tg^p- This means that for a star 
cluster to reach a significant fraction of the Hill radius of 
the nuclear ring is very fast, i.e., typically of order less 
than 1% of the dynamical friction time. This implies a 
factor 20 — 30 quicker than tgcp, or, with icross = 20 Myr, 
less than 20 Myr for q' = 0.01 and less than 2 Myr for 
q' = 0.1. We come back to this separation in the next 
section when we (numerically) investigate the system in 
greater detail. 

3.2. Separation and Density Enhancement 

We solve equations pi)) and (fTSl) using standa rd ex- 
plicit finite-difference methods (jPress et al.l 11993 ). For 
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Fig. 2. — The top panel shows the evolution of the radial position 
of a star cluster, rs, relative to the initial outer radius of the nuclear 
ring, r^j Q, as function of the time, t, in units of the initial dynam- 
ical friction time, 4df,o- The bottom panel shows snapshots of the 
surface density profile at the edge of the nuclear ring for various 
times, separated by At/tdf,o = 0.06, with labels A-D as indicated 
in the top panel. The star cluster has mass of 1% (g' = 0.01) 
relative to the gas mass in the nuclear ring. Furthermore, the di- 
mensionless model parameters are set to a' = 0.01, /9' = 1, and 
7' = (no mass inflow). 



good spatial resolution, we use 2000 grid points in r/r^.o 
between zero and 1.50, with the star cluster initially at 
fs/fdfi = l-Ol- We set Qd = Md/Mcnc — 0.01 similar 
to the observed ratios of (gas) mass in the ring to en- 
closed mass within the ring radius derived in § 12.11 Fur- 
thermore, for the dimensionless parameters a' and j3' , 
representing the strength of the ring viscosity and tides 
relative to the dynamical friction, we choose /3' = 1 and 
a' = 0.01. We first study the case without source term, 
i.e., 7' — 0, while in the next section § 13.31 we include 
the effect of gas mass inflow. 

In Figure [21 we show the (numerical) solution of equa- 
tions (fT^ and ([T5)) for a star cluster with mass relative 
to the (gas) mass in the nuclear ring of q' — 0.01. The 
top panel shows the evolution of the radial position of 
the star cluster, r^, relative to the initial outer radius of 
the nuclear ring, r^fl, as function of the time, t, in units 
of the initial dynamical friction time, idf,o- The bottom 
panel shows snapshots of the surface density profile at 
the edge of the nuclear ring for various times, separated 
by Ai/idf.o = 0.06, with labels A-D as indicated in the 
top panel. The initial built-up of material at the edge of 
the ring is very fast and only after i/idf.o — 10~^ (point 
A) an enhancement of about 30% in the surface mass 
density is achieved just inside Tdfl. This rapid initial en- 
hancement is due to the fourth-power radial dependence 
of the tidal torque in eq. (|13p . and may not develop in 
a realistic system. Nevertheless, the initial ring tidal 
torque will be very strong and pushes the star cluster 
outward beyond the initial edge of the ring. After a time 
i/idi,o = 0-06 (point B), the star cluster has moved out 
about 20% from its initial position. The surface mass 
density near the outer edge decreases now with respect 
to its initial value due to viscous spreading. Points C and 
D continue the trend of an outward moving star cluster 



Fig. 3. — Same as Figure [2] but for at ten times more massive 
star cluster with q' = 0.1. After an initial outward migration 
the star clusters turns around and starts moving inwards. At the 
same time, the surface mass density is pushed to smaller radii and 
enhanced by a factor three to two. 



and viscously spreading ring, albeit at a decreasing pace. 

In contrast, for the case of g' = 0.1 (a ten times more 
massive star cluster), the star cluster does not contin- 
uously move outward, but rather stop and even moves 
inward as shown in Figure [3] As in the earlier case, 
an initial surface mass density enhancement is built at 
point A, though it is a factor of ten larger than for the 
case q' = 0.01, i.e., 300% {q' = 0.1) as opposed to 30% 
{q' = 0.01) of its initial value. However, the subsequent 
evolution of this case differs from the earlier case. The 
star cluster still moves outward, but it reaches only about 
15% of its initial position. Its maximum in rs/rd,o is 
reached between point B and C, after which the star 
cluster "turns around" and starts moving inward with 
respect to the initial outer radius of the ring. The sur- 
face mass density at the edge of the ring monotonically 
moves inward, effectively being "shepherd" inward by the 
star cluster, similar to the scenario in COS. 

In both cases, the radial separation between the star 
cluster and the (edge of the) nuclear ring approach a con- 
stant value, which we expect from the order of magnitude 
estimate in §[2T2]to be close to the Hill radius of the ring, 
TH.d, as defined in eq. ^. Substituting qd = 0.01 and 
using that r^ > Vdfi ^ i"H,d, we estimate this Hill ra- 
dius to be around one-fifth of the initial radius of the 
ring, i.e., rH,d/T'd.o ^ 0.2. Combining the above shifts in 
the radial positions of the star cluster and nuclear ring, 
we find that the separations for both cases indeed con- 
verge to around one-fifth of rd,o, although for q' = 0.01 
the separation stays slightly smaller, while for q' = 0.1 
it is somewhat larger. This is further illustrated in Fig- 
ure |4l where we show the evolution of the radial separa- 
tion rs — Td.edge normalized by the Hill radius rH,d- We 
use two definitions of the radius of the disk edge, r dredge- 
For q' — 0.01, we define rd^edge as largest r^ for which 
S/So > 75%, while for q' = 0.1, rd,edge is defined as 
the rd for which S peaks. Despite these two different 
definitions. Figure 0] shows two important points. First, 
the separation between the disk and satellite is of order 



Migration of Star Clusters and Nuclear Rings 



0.5 




0.2 



Fig. 4. — The evolution of the radial separation rs — Td^^dge 
between the star cluster and the edge of the nuclear ring, normal- 
ized by the Hill radius r^ ^ of the nuclear ring. The top curve is 
for a ten times more massive star cluster (q' = 0.1, see Figure [Sj 
then bottom curve {q' = 0.01, see Figure [Jjl. For both cases, the 
time to converge to the final separation is about 20% of the initial 
dynamical friction time t^f q, while the initial separation is very 
fast, reaching half of the final separation in less than 1% of tdf,o- 



fH4-, consistent with our expectations from § 12.21 Sec- 
ond, the final separation is reached after about 20% of 
the dynamical friction time, consistent with the estimate 
of isepAdf,o ^ TH.d/rdfi ^ 0.2 in § [HH] for qd = 0.01. 
This figure also confirms that the initial separation is 
very quick, indeed reaching half of the final separation 
within a one percent of the dynamical friction time, also 
consistent with our earlier expectations. 

3.3. Constant Mass Inflow Rate 

We now include the effect of a constant gas mass infiow 
on the dynamics of the cluster-ring systems. Figure [5] 
shows the results of (numerically) solving the same sys- 
tem with g' = 0.1 as in Figure [3l except that the source 
term in eq. p4p now contributes with 7' = 1. We es- 
timated in § 13.11 that the latter corresponds to an mass 
inflow rate of M ^ O.lMQyr^^. Comparing the top 
panels of the two figures, we find that the star cluster 
moves (slightly) further out when 7' = 1 and remains 
constant at the maximum radial position, as opposed to 
the shrinking of r^ after having moved out for an initial 
period when 7' = 0. This is because in the case that 
7' = 1, the inflow torque is balancing the torque due 
to dynamical friction, as we estimated in eq. ([7]). The 
bottom panels of the two figures are very similar, both 
showing an enhancement as well as inward migration of 
the surface mass density. The most significant difference 
is that in Figure[5]beyond the edge at r/rdfi ~ 0.92 there 
is a long tail out to r/r^fi ~ 1, due to mass inflow into 
this region. Another difference is that for 7' = 1 after 
the initial inward push of the edge, it remains flxed at 
T/fd.o ^ 0.92, whereas for 7' = it still moves inward, 
although at a decreasing rate. This is consistent with the 
above difference in the change in the radial position r^ of 
the star cluster. As a consequence, also in the case that 
7' = 1, the separation between the star cluster and the 
(edge of the) nuclear ring converge to a constant value. 
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Fig. 5. — Same as Figure |3] with q' = 0.1, but now with 7' = 1. 
The torque due to the latter mass inflow rate approximately cancels 
the dynamical friction torque or the star cluster, preventing its 
radius r^ to shrink. The long tail extending from r/r^ q ^ 0.92 
out to r/rd,o '^ 1 is due to mass inflow into this region. 
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Fig. 6. — The evolution of the radial position of a star cluster 
with q' = 0.1, relative to the initial outer edge of the ring, r^ — 
r^j Oi in units of the Hill radius of the ring, r^ ^j, for four different 
constant mass inflow rates 7' = 0, 0.3, 1 and 3 from bottom to 
top. After a similar initial outward migration, the motion of the 
star cluster goes from in-falling for 7' < 1 to outgoing for 7' -^ 1. 



which is again of the order of the Hill radius of the ring. 
In Figureini we show the evolution of the radial position 
of the star cluster relative to the initial outer edge of the 
ring Ts — rd.Oj in units of the latter Hill radius rH,d, for 
four different constant mass inffow rates. From bottom 
to top the curves correspond to 7' = 0, 0.3, 1 and 3. 
This clearly shows the transition around 7' = 1: after a 
similar initial outward migration, the motion of the star 
cluster goes from in-falling for 7' < 1 to outgoing for 
7' > 1. This is fully consistent with our estimated scale 
for the mass infiow rate M in eq. ([7]), case of single star 
cluster with q' = 0.1. 
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4. COMPARISON TO OBSERVATIONS 

We describe the expected evolution of nuclear rings due 
to interaction with star clusters based on our findings 
in the previous sections. First, we focus on the evolu- 
tion in radius, and discuss possible azimuthal variations 
thereafter. Finally, we take a more detailed look at the 
well-studied nuclear ring in NGC 4314. 

4.1. Radial Evolution of Nuclear Rings 

From our analysis in §[51 we expect that a newly formed 
star cluster moves outward until it is about a Hill radius 
separated from (the edge of) the nuclear ring. We de- 
rived in ii l2.2l a range in rn^d from about 50 pc to 500 pc, 
with 150 pc as a typical value. For a host galaxy at 
distance of ^ lOMpc, this corresponds to about 1" to 
10", and typically 3", which should be readily observable. 
However, while the position of the (young) star clusters 
with respect to the center of their host galaxy is straight- 
forward to measure, the radial extent of the gas in the 
ring is often difficult to determine. It requires high spa- 
tial resolution CO measurements to trace the underlying 
molecular gas. Also, part of the gas (at the edge) of the 
ring may be severely disturbed due to winds that clear 
out the gas and dust from the newly formed star clusters 
and/or ionized by the massive, bright stars in the young 
star cluster (but see ^ 15. 31 and ? l5.4p Still, the dust seeni s 
to be well mixed into the gas (e.g. iBrunner et aI1l2008D . 
indicating that dust provides a good tracer of the extent 
of the gas in the ring. Apart from the contact points 
where the bar dust lanes connect with the nuclear ring, 
the dust is indeed found on the inside of the distribution 
of the star clusters (see also § 12. 2|) . 

Still, using the (peak in the) dust distribution as an es- 
timate of the edge of the ring, the (unique) measurement 
of a separation is complicated by the spread in radius of 
the star clusters. Part of this spread is because of the 
dependence of the dynamical friction time on the mass 
of the star cluster: tdf/^cross — l/<z'- Even though we 
found that it takes less than 1% of the dynamical fric- 
tion time to reach about half of the Hill radius, with 
a typical crossing time of tcross = 20 Myr, this means 
only 2 Myr for the most massive star clusters {q' = 0.1), 
but as much as 200 Myr for the least massive clusters 
(g' = 0.001). Since i^opAdf ^ 0.2 - 0.3, to reach the 
full separation takes a factor 20 — 30 longer. This means 
that only the most massive clusters reach at least half 
of the Hill radius separation before all bright, massive 
stars die. The less massive clusters have already faded 
significantly, or even dissolved. So, we then might expect 
a smooth radial gradient with the most massive clusters 
on the outside, which in turn can be used to establish 
the Hill radius. 

Unfortunately, we found that the more massive cluster 
indeed rapidly reach a maximum separation, but after 
staying a while at this radius (depending on the mass 
inflow rate), they move inward again by pushing the 
ring inward (compare for example the top panels of Fig- 
ures [2] and [3]) . Hence, over time the less massive clus- 
ters, if they survive, may catch up with the more mas- 
sive clusters. We conclude that only among the (very) 
young (few Myr) clusters we might expect a radial gra- 
dient with the more massive clusters to the outside, from 
which the Hill radius can be estimated. Later on, if a ra- 



dial gradient is visible at all, it might even be opposite, 
with the more massive clusters on the inside. Since the 
more massive clusters have higher survival chance, the 
latter 'turn around' may result in at most a mild radial 
gradient i n age (wit h the older clusters on the inside). 
JMazzuca et al.l (|2008i) find only for 2 (out of 22) nuclear 
ring galaxies, a radial age gradient with indeed (slightly) 
older clusters on the inside. 



4.2. Azimuthal Age Gradient 
Unlike the radial age gradient. iMazzuca et al.l 



as in other detailed studies (e.g. by lB5ker et ahl 



2008 ) 



2om 



do find signatures of an azimuthal age gradient in most 
of the nuclear rings, ranging from clean bipolar age gra- 
dients to a correlation between the youngest star clusters 
and the contact points of the bar dust arms with the nu- 
clear ring. The natural explanation is that at (or very 
close to) these contact points the inflowing gas shocks 
and forms stars, after which the resulting star cluster 
travels in the ring away from the contact points. Indeed, 
for a typical b ar pattern speed of ~ 50 km s^ ^ kpc~ ^ (e.g. 
iGerssenI l2002f) the corresponding time for the contacts 
points to rotate is ^ 120 Myr, significantly shorter than 
the local crossing time icross ^ 20 Myr. Alternatively, 
if the star formation occurs due to (local) gravitational 
instability of the gas in the ring, the star clusters are ex- 
pected to form throughout the ring without an azimuthal 
age gradient. Instead of associating th is instability wit h 
an inner Lindblad resonance (ILR; e.g. lElmegreenlll994f ). 
it is a natural consequence of shepherding: The gas in 
the ring loses angular momentum to the star cluster(s), 
which in turn lose angular momentum to the host galaxy 
through dynamical friction. This means that gas at the 
edge of the ring moves inward, resulting in an increase in 
the gas surface density at the edge (see also Figures [51 [21 
and [5]), until its becomes unstable against self-gravity. 

In the case of shepherding we do expect, as observed, 
a mixture in azimuthal age gradients. Initially, when 
the gas surface density in the ring is below the criti- 
cal value for star formation, star clusters only form in 
episodes at the two contact points, they move together 
with the remaining gas along the ring, resulting in bipo- 
lar age gradient. The newly formed star clusters move 
outward, while at the same time 'pushing' the gas in- 
ward, causing an over-density built-up at the edge of the 
ring. When the gas surface density somewhere along the 
edge of the ring reaches the critical value, it becomes un- 
stable and forms stars. The more massive the star cluster 
(or group of star clusters acting together), the higher the 
over-density and the sooner the critical value is reached. 
Generally, the most massive clusters and hence the star 
formation might be anywhere along the ring. However, 
depending on the mass inflow rate, there might be rela- 
tively more gas near the contact points to start with, so 
that first instability-driven star clusters still form near 
the contact points. Therefore, an azimuthal age gradient 
might still exist over part of the ring, but eventually the 
gradient is expected to flatten. 

We found in § 13.31 that for a low enough mass inflow 
rate, the (most massive of the) existing star clusters move 
inward by dynamical friction, while 'pushing' and in- 
creasing the surface mass density of the gas in the outer 
part of the nuclear ring. As a result, the gas might reach 
the critical surface mass density and form new stars (in 



Migration of Star Clusters and Nuclear Rings 



clusters) all along the edge of the ring. At the same 
time, near the contact points, which also move inward 
along with the nuclear ring, star formation from inflow- 
ing and shocking gas might be ongoing. We argue below 
in ii l4.3l that the nuclear ring in NGC 4314 is an example 
of this mechanism. Since even in the case of a high mass 
inflow rate, eventually enough (massive) star clusters are 
formed for their (combined) dynamical friction to 'push' 
to ring inward, the corresponding gas and bright star 
clusters migrate inward from larger to smaller X2 orbits. 
In other words, nuclear rings will be located between the 
so-called outer and inner ILR, although strictly speaking 
resonances are only applicable to systems with a weak 
bar. Depending on the details of the mass inflow rate 
and star formation efficiency over time, we thus expect 
a variety in sizes of nuclear rings with respect to the res- 
onance radii, as well as a significant diversity in their 
star cluster distribution. Hydrodynamic simulations in- 
deed show the inward migration of nuclea r rings (e.g. 
iReganfc TeubenI l2003t iFukudaet aLll200 0V and obser- 
vations of nuclear rings c onfirm the diversity i n their star 
cluster distribution (e.g. iMazzuca et al.|[2008[ ) 

4.3. Nuclear Ring in NGC 4314 

NGC 4314 is a nearby {D - lOMpc, 1" ~ 48.5 pc) 
galaxy with a large-scale bar out to a radius of ~ 65" 
(3.15kpc). It host a prominent circumnuclear ring of 
star formation that peaks at a radius of ~ 7.0" (340 pc), 
that is visible in Ha, radio continuum as well as optical 
color maps. Hi g h spa tial resolution CO observations by 
[Benedict et al.l ()1996f ) clearly show that the molecular 
gas peaks at a smaller radius of ~ 5.5" (270 pc). 

The circular ve locity curve flattens at ^ 7" to a value 
of ~ 190kms~^ (jCombes et al.l [199^ . corresponding to 
a (spherically) enclosed mass Menc ~ 2.8 x 10^ M©. The 
total molecular gas mass within 7" is ~ 15 x IO'^Mq. 
However, this includes gas that is still flowing in along 
the dust lanes at the contact points, as well as a signif- 
icant amount in the center where a strong peak in the 
CO coincides with a peak in Ha + [Nil] consistent with 
the LINER definition. Hence, we estimate the r ing-to- 
enclo sed mass to be qd ^^ 0.03 (see also Benedict et al.l 
Il996t ). The Hill radius of the ring is thus rn ^^ 105 pc or 
rn ~ 2.2". The observed difference of ~ 1.5" between 
the peaks in star formation and CO flux is smaller than 
the predicted Hill radius, which is expected based on the 
discussion in § 14.11 Although the peak in CO flux nicely 
corresponds with the peak in the dust distribution and 
thus likely traces the edge of the ring, the star clusters 
have not yet reached the Hill radius. Given a separation 

1 /3 

time of tsep/idf ^ Qd ^ 0.3, a dynamical friction time 
of td{/t cross — 1/g', a crossing time of Across '^ llMyr, 
and an average age of the star clusters of i ~ 15 Myr, 
we expect based on eq. (|T8l) . that the star clusters have 
reached a separation of about Ar ~ l.bq'^^^rn.d- For 
star clusters with a (combined) mass relative to the gas 
mass in the ring from q' = 0.01 to q' = 0.1, this means 
a separation from about 1.0" to 1.8", bracketing the ob- 
served difference. 

There is no evidence for an azimuthal age gradient, but 
a slight tendency for a radial age gradient with the few 
clusters outside the ring begin all older than 15 Myr. This 
is consistent with the shepherding scenario in which the 



gas in the ring has been 'pushed' inward until a critical 
surface density is reached and all along the edge of the nu- 
clear ring star clusters are formed, which then move out- 
ward. The (group of) star cluster(s) responsible for the 
inward migration of the ring likely have faded too much 
to be visible against the old underlying stellar popula- 
tion of the host galaxy. However, further support for the 
inward migration of the ring comes from two blue stellar 
spiral arms that seems to connect the outer-most X2 or- 
bits ( also referred to as th e outer IL R) with the nuclear 
ring (jWozniak et al.l Il995t [Benedict et all |2002|). Sub- 
tractin g a model of the red host galaxy. [Benedict et al.l 
(|2002l ) find that the stars in the arm close to the nu- 
clear ring are ~ 20 Myr old, with an indication of an 
increase in age beyond ~ 100 Myr going outward in the 
arms. Whereas during the inward migration the gas sur- 
face density along (the edge of) the ring was insufficient 
to form stars, shocks at the contact points might still 
have allowed the formation of stars. While the contact 
points moved inward together with the ring, these stars 
st ayed behind an d creat ed the blue arms. 

[Benedict et al.l ([1996[ ) derive an outer ILR at ~ 13" 
and inner ILR at ~ 4" for an assumed bar pattern speed 
Vtp ^ 72kms~^ kpc^^ of the large-scale bar. The latter is 
likely an upper limit because typically the co-rotation ra- 
dius is beyond the ~ 3 kpc radius of the bar, so that with 
a 190kms~^ circular velocity, fip < 60kms^^ kpc^^. In 
turn this implies that the above outer and inner ILR 
are probably lower and upper limits respectively. This 
clearly shows that it is critical to obtain an independent 
measurement of the bar patter n speed, using for example 
the Tremaine- Weinberg ([19841 ) method. Even so, the nu- 
clear ring is in between the outer and inner ILR, in line 
with the shepherding scenario, in which the ring migrates 
inward until the increasing gas surface density at the edge 
reaches the critical value for star formation. While the 
newly formed star clusters move (slightly) outward peak- 
ing at r^ 7.0", the gas in the ring is 'pushed' further in - 
ward peaking at ~ 5.5". Note that [Benedict et al] ([1993[ ) 
detected a nuclear bar within ~ 4" that might provide 
a mechanism to transport part of the gas further inward 
towards the presume d black hole in the center (see also 
[Benedict et"al[[2002[ ). 

5. DISCUSSION 

We discuss some of the assumptions of our model and 
we speculate about possible consequences the inward mi- 
gration of the cluster-ring system has on transporting gas 
to the center of the galaxy. 

5.1. Isothermal Density Profile 

In deriving the equations for our model in § 12. 3( 
we have adopted a density profile p oc r~^ for the 
(bulge of the) spiral galaxy in which the star cluster 
and nuclear ring are embedded. At the radii wher e nu- 
clear rings are present (0.2-1.7kpc: [Mazzuca et al"l2 008). 
the density of galaxies indeed seem to be not too far 
from such an isothermal profile, as indicated, for ex- 
ample, by studies that combine stellar dynamics with 
strong gravitational len sing (e.g. [Koopmans et al.l[2006l : 
[van de Ven et al.|[2008[ ). A more straightforward way is 
to consider that the surface brightness profiles of (bulges 
of) galaxies are well fitted by a Sersic (1968) profile 
I{R) (X exp[— (i?/i?e)^/"], with index n and the effective 
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radius Re enclosing half of the total light. The (numer- 
ically) de projection is well approxi mated by the density 
profile of iPrugniel &: SimienI ()1997f ). of which the (nega- 
tive) slope is given by 



n \Re 



l/n 



(19) 



with p„ = 1.0 - 0.6097/n + 0.05463/n^ and &„ = 2n - 
1/3-H 4/405 n-H 46/25515 n^ (ICiotti fc Bertinl[l999l) . For 
an a verage bulge of a spiral galaxy n ~ 2 and i?e ~ 1 kpc 
(e.g. iHunt et al.ll2004D . so that at the typical radius r ^ 
0.5 kpc of a nuclear ring, we indeed find that the slope of 
the density 7 ^ 2 is close to isothermal. 

5.2. Gas Clumps and Clouds 

The star cluster that is being formed in the nuclear 
ring is already subject to dynamical friction while it is 
still a gas clump, as long as the material that makes up 
the (proto) star cluster is bound enough relative to the 
background stellar population of the host galaxy. Since 
this gas clump contains (significant) more mass than the 
final cluster of only stars, the initial dynamical friction 
is (much) stronger. However, the clearing of gas (and 
dust) from the newly formed star cluster is very fast and 
efficient, since even the yo ungest visible star clusters are 
only mildly reddened fe.g. iMaoz et al.|[2001h . 

We assume that the gas in the nuclear ring is uniform, 
i.e., no gas clouds (or clumps) aside from those which 
turn into star clusters. Since such gas clouds would also 
be subject to dynamical friction from the background 
bulge stars and thus move inward, dynamical friction on 
the nuclear ring itself has also been proposed as an ex- 
planation of the inward migration of the nu clear ring in 
for example NGC4314 (|Combes et al.|[T993) . 

Suppose that the gas is made up of N gas clouds of 
the same mass Mc = M^/N, then about N = 1/q'^ 
clouds are needed to have a combined dynamical fric- 
tion torque of the same order as that of a single star 
cluster of mass Ms = q'Md- Adopting as before the 
two cases of q' = 0.01 and q' — 0.1, this respectively 
means N = 10* and N ~ 10^ gas clouds of mass 
Mc = 5 X 10^ Mq and M^ = 5 x 10^ Mq for a typical nu- 
clear ring mass of Md = 5 x 10^ Mq. From § l2.1l we have 
that the dynamical friction time in units of the cross- 
ing time is approximately given by idf /Across — Md/Mg, 
so that with Ms — Mc = Md/N for a gas cloud we 
obtain idf — Ntcmss- Given a typical crossing time 
of Across ^ 20 Myr, the timescales of dynamical friction 
for the latter two cases are about tdf ~ 200 Gyr and 
idf ^ 2 Gyr, respectively. We thus conclude either nearly 
no inward migration [q' ^ 0.01), or gas that is very 
clumpy with quite large gas clouds {q' ~ 0.1). 

Moreover, since the gas clouds will span a range in 
mass with different corresponding dynamical friction 
timescales, they migrate inward at different rates, leading 
to a (radial) diffusion of the nuclear ring, especially with 
respect to the gas that is not in clouds and stays behind. 
Although in some galaxies, like M 83, that are claimed to 
harbor a (partial) nuclear ring, the corresponding mor- 
phology is rather clumpy and chaotic, most of the nuclear 
rings are are quite smooth and well defined, including the 
nuclear ring in NGC4314. CoUisionless star clusters are 
also less prone to destruction than coUisional gas clouds. 



in particular in the presence of shear stresses from dif- 
ferential speeds. As a result, at least the most massive 
star clusters will sink faster inward than the gas clouds, 
and appear to be on the inside of the gas (clouds) in the 
nuclear ring, which typically is not observed. 

5.3. Star Formation within the Nuclear Ring 

We have assumed that the star cluster form just exte- 
rior of the nuclear ring. This assumption simplifies the 
physics and the numerical computation, though it is not 
absolutely necessary. If the star cluster forms near the 
outer edge of the nuclear ring, we do not expect the re- 
sulting dynamics to be significantly different. By near, 
we demand that the separation between the star cluster 
and the outer edge is within one Hill radius of the nuclear 
ring. In this way, the amount of mass which participates 
in the tidal evolution of the cluster-ring system exterior 
to the star cluster is less then interior to the star cluster. 
As long as this condition is fulfilled, the evolution should 
proceed along the lines of our calculation. 

There are a few reasons to believe that the star clus- 
ters (preferably) form near the outer edge of the nuclear 
ring. The supply of gas from the bar merges with the 
ring at the outermost radius. It is unknown how the gas 
spreads radially over the rest of the ring, but if it is due 
to viscous processe s, for insta nce the magneto radial in- 
stability (MRI; Bal bus &: Haw lcv 1991), then the surface 
density is maximal at the outer edge. As a result, the 
gas in the ring is most likely to be(come) unstable near 
the outer edge and hence preferentially fragment there. 

While many clusters will form near the outer edge 
(and hence susceptible to shepherding) for the reasons 
given above, clusters are also likely to form inside of the 
outer edge. For instance, once massive star clusters are 
formed, their interaction with the gas ring may result 
in large and fairly symmetric enhancements in the gas 
surface mass density (see for example Figure [3]), so that 
new, but significant smaller, star clusters can be formed 
around these enhancements. For the relatively small star 
clusters formed on the inside, we expect the inward mi- 
gration to follow that of forming proto planets (see e.g. 
IWard 1997). In addition, if the gas flowing along the bar 
is very clumpy, possess a different vertical scale height, 
or a different inclination than the gas ring, it may not 
(only) merge at the outermost radius. In this way, even 
rather large star clusters might form on the inner edge 
of the nuclear ring. They then migrate further inward 
(relatively faster than cluster formed on the outer edge 
as they encounter less of the gas ring), possibly along an 
inward extension of the inspiraling gas. Whereas these 
are possible explanations for the observation that some 
star clusters appear more on the inside of the nuclear ring 
(see also § 12. 2[) , alternatively they might form already in- 
side of the ring as the result of the local interaction of 
a non-axisymmetric structure such as a nuclear bar or 
spiral with the interstellar medium (see also § 15. 5[) . 

5.4. Stellar Feedback 

Winds originating from (massive) stars and/or super- 
novae in the star clusters might not only clear out the 
remaining gas in the star cluster itself, but also push 
away part of the gas in the nuclear ring. This creates a 
separation between the cluster and the ring that might 
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mimic the outward migration of the star cluster. How- 
ever, since the density of the interstellar medium is far 
less than that of the gas in the ring, the winds will es- 
cape in the vertical direction like a fountain arising from 
the (equatorial) plane of the ring. Therefore, even if the 
cluster is embedded in the ring itself, the winds are ex- 
pected to clear out at most a region that is equal to scale 
height of the disk, hd, which is significantly smaller than 
the Hill radius of the ring, r^ii- 

To show this, we start from iToomrd ()1964[ ) Q pa- 
rameter for the gas in the ring: Q = Csk/ttCS 
(jBinnev fc Tremaiii3ll987[ ). Here, the epicycle frequency 
K is defined as k^ = {4 + dlnfi/dlnr)!]^, so that k = 
\/2^ for the assumed (isothermal) density profile p oc 
r^^. Adopting Cg — hj^, we obtain for the scale height of 
the ring hd/rd — (Md/Mcnc) Q/V^- On the other hand, 
the Hill radius of the ring rH,d/rd ^ (Md/Me„c)^/^. As- 
suming Q — 1 and substituting the mean observed ring- 
to-enclosed mass ratio Md/Mcnc — 0.02, we find that 
hd/i'H,d — 0.05. This indeed shows that winds might 
clear out a region that is at most a few per cent of the 
Hill radius of the nuclear ring. 

A second possible stellar feedback effect, is due the (ul- 
traviolet) radiation coming from the massive, bright stars 
in the (young) star cluster, which might (photo)ionize the 
(molecular) gas in the nuclear ring. We can estimate the 
radial extent of this ionization from the radius rion of the 
Stromgren sphere. We start from iVjon — ari^iTrrf^^/S, 
with cross section a ~ 3 x 10~^^ cms~^. We assume that 
bright stars radiate roughly at 10% of the Eddington lu- 
minosity in the ultraviolet (UV; > 10 eV), to obtain an 
ionization rate of TV < 10*^ (M^^/Mq) photons per sec- 
ond. With the above estimate of the scale height hd, 
we find for the gas mass density pg ~ Afonc v2/7rrj^(5. 
Substituting Afcnc ~ 2.5 x 10^ Mq, r^ ~ 0.5 kpc and 
Q ~ 1, and dividing by the proton mass, we find for 
the number density of the gas Ug ^ A x lO^cm"'^. In 
this way, we estimate that a (UV) radiating star of mass 
M^ is able to ionize gas out to a (maximum) radius of 
no„<3(A/,/M0)i/3pc. 

For a typical O star of M^, ~ 40 M© , this means non < 
10 pc, which is significantly smaller than the typical Hill 
radius of rn^d ^ 150 pc. The combined (and focused) 
ionizing radiation of thousands of massive young stars 
with a lifetime of only a few to ten Myr are needed to 
reach the latter radius. Even then, the dust, which is 
hardly effected by the radiation, will still provide a good 
tracer of the gas in the ring, and as such can be used 
to establish the separation between star clusters and the 
nuclear ring. 

5.5. Feeding the Nucleus? 

For sufficiently massive star clusters and/or for a suf- 
ficiently large number of them, it is in principle possible 
for these star clusters to shepherd the nuclear ring to 
smaller radii. As in COS, we assumed the gas ring sits 
in a simple spherically symmetric potential. However, 
the gas ring arises as a result of the stellar bar, which 
is a non-axisymmetric feature in the potential, driving 
the gas inward towards the Xi — X2 orbit transition. In 
additi on, the gas ring is self-gravitating with Toomre's 
(|1964[ ) Q near unity. We now speculate on the effect of 



herding scenario. 

First, as the gas ring is being shepherded inward, it 
moves from one X2 orbit to a smaller one. Eventually, 
it may move inward of the sma llest non-intersecting X 2 
orbit in a barred potential (e.g. iRegan fc Teubenl l2003'). 
The gas which is pushed beyond this point will find it- 
self in intersecting orbits, losing angular momentum on 
each orbit as it shocks against gas on the same orbit or 
on neighboring orbits. The high(er) density regions re- 
sulting from the shocks, might give rise to the observed 
(chaotic) spurs of gas and dust in the centers of (late- 
type) galaxies (e.g. lElmegreen et al.lll998l : lMartini et al.l 
[2003fl . 

Second, the effect of self-gravity in the shepherded gas 
ring is nontrivial. As we stated above in § 15. 3i the ef- 
fect of an enhanced surface density at the edge of the 
ring results in Q < 1, if the initial Q ~ 1. In addi- 
tion, this local enhancement in the surface density may 
give rise to the gaseous analogy of s tellar self-gravitating 
groov e and edge modes (see a lso Lovelace & Hohlfcld 
1978t ISellwood fc KahnI il991l: Papaloizou fc Savonije 
This may result in the formation of nu- 



1992). 



clear bars (e.g. iLaine et al.ll2002t lE rwin fc Sparki 120021: 
Englmaier fc Shlosman 12004 ) or nuclear spira l s (e.g . 
Englmaier fc Shlosmanl I2OO0I: iPogge fc Martinil l200a) . 



this non-axisymmetry and self-gravity on the gas shep- 



providing a means for inward transport of gas from the 
nuclear ri ng region in analogy to the "bars within bars" 
scenario (jShlosman et al.|[l990[ ). 

6. CONCLUSIONS 

We have studied the dynamics between a star cluster 
and the gas in the nuclear ring from which it was formed. 
The star cluster is subject to dynamical friction with the 
background bulge stars in the host galaxy. The resulting 
dynamical friction torque induces a tidal (and viscous) 
torque on the gas in the nuclear ring. In addition, we 
investigated a torque due to in-falling gas along the bar 
onto the nuclear ring at the contact points. We arrived at 
the following, observationally testable predictions based 
on our model. 

(i) The star clusters which are formed near the outer 
edge of the ring migrate outward. The final sepa- 
ration is of the order of the Hill radius of the ring, 
which is 20 — 30% of its radius if the ring's gas mass 
is the observed 1 — 3% of the enclosed mass. Simi- 
larly, the time to reach this final separation is about 
20 — 30% of the dynamical friction time. The initial 
separation is very fast, reaching half the Hill radius 
within only 1% of the dynamical friction time. The 
latter can be as short as a few million years for a 
massive enough (> 10® Mq) star cluster. 

(ii) While such a massive cluster moves out the surface 
mass density of the gas at the edge of the nuclear 
ring is pushed inward and gets enhanced by a fac- 
tor of a few. If this enhancement reaches the crit- 
ical value for star formation, new star clusters are 
expected to form along the ring. This is in addi- 
tion to the (ongoing) formation of star cluster from 
shocking gas near the contact points. The (initial) 
azimuthal age gradient that is expected from the 
latter gets diluted or even washed out by the for- 
mer (secondary) star cluster formation. 
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(iii) If the (formed) star cluster are massive and/or nu- 
merous enough, and the gas mass inflow rate low 
enough (or even shut down), we expect the inward 
migration of the cluster-ring system as a whole. 
This mean that, even though initially the nuclear 
ring could have formed as often suggested at an 
(inner) Lindblad resonance, it can migrate closer 
to the center of the host galaxy. 

The latter implies that the most massive star cluster, 
after an initial outward separation, turn around and mi- 
grate inward along with the nuclear ring. Consequently, 
any initial radial age gradient due to the dynamical fric- 
tion time being inverse proportional to the mass of the 
star cluster, is diluted or even inverted. Even so, both a 
radial and azimuthal age gradient are observationally dif- 
ficult to measure. When the most massive stars are still 
alive, the bright emission lines can be used for (relative) 
age dating (e.g. iBoker et al.l[2008( l. but after a few tens 
of million years it becomes already very challenging to 
distinguish the fading star clus ters from the backgro und 
stellar population (but see e.g. [Benedict et aLll2002f ). 

The outward migration of a least the most massive 
star clusters is easier to establish, and using the dust to 
trace the gas, observations indeed show the young star 
clusters exterior to the nuclear ring. The measured sep- 



aration should be of the order of the Hill radius of the 
nuclear ring, which indeed is the case in NGC 4314. This 
well-studied nuclear ring also seems to have migrated in- 
ward, leaving a tail of older star clusters behind and be- 
ing currently located between the inner and outer ILR. 
To firmly confirm such a "shepherding of gas" towards 
smaller radii, the dynamics of the host galaxy, including 
the (bar) pattern speed, needs to be measured accurately. 
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